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Abstract 

cn 

We develop two inverse scattering schemes for locating multiple electromagnetic 
(EM) scatterers by the electric far-field measurement corresponding to a single in- 
cident/detecting plane wave. The first scheme is for locating scatterers of small size 
compared to the wavelength of the detecting plane wave. The multiple scatterers 

I ^ could be extremely general with an unknown number of components, and each scat- 

terer component could be cither an impenetrable perfectly conducting obstacle or a 
penetrable inhomogeneous medium with an unknown content. The second scheme 
is for locating multiple perfectly conducting obstacles of regular size compared to 
the detecting EM wavelength. The number of the obstacle components is not re- 
quired to be known in advance, but the shape of each component must be from a 
certain known admissible class. The admissible class may consist of multiple differ- 
ent reference obstacles. The second scheme could also be extended to include the 
medium components if a certain generic condition is satisfied. Both schemes are 
based on some novel indicator functions whose indicating behaviors could be used 

i i to locate the scatterers. No inversion will be involved in calculating the indicator 

functions, and the proposed methods are every efficient and robust to noise. Rigor- 
ous mathematical justifications are provided and extensive numerical experiments 
are conducted to illustrate the effectiveness of the imaging schemes. 

1^ 1 Introduction 

This paper is concerned with the inverse scattering problem of reconstructing an in- 
homogeneous scatterer located in an otherwise homogeneous space by measuring the 
corresponding electromagnetic (EM) wave fields far away from the scatterer produced 
by sending some detecting EM wave fields. The inverse electromagnetic scattering prob- 
^ lem has been playing a key role in many areas of science and technology, such as radar 
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and sonar, non-destructive testing, remote sensing, geophysical exploration and medical 
imaging to name just a few; see (4[[5j[Io|[TT}[T5j[20j[25}[3TJ[33] and the references therein. 
In the current article, we shall mainly consider the reconstruction scheme for this in- 
verse scattering problem. There are extensive studies in the literature in this aspect and 
many imaging schemes have been developed by various authors, and we would like to 
refer to [T}^[^[^[g[l7}{l9}[2T}[30}[32] and the references therein. However, we would 
like to remark that in those schemes, one either needs to make use of many wave mea- 
surements or if only a few wave measurements are utilized, then one must require that 
the underlying scatterer is of small size compared to the detecting wavelength. In this 
work, we shall consider our study in a very practical setting by making use of a single 
electric far- field measurement. That is, we shall consider the reconstruction by measur- 
ing the far- field electric wave corresponding to a single time-harmonic plane wave. From 
a practical viewpoint, the inverse scattering method with a single far-field measurement 
would be of significant interests, but is highly challenging with very limited progress 



in the literature; we refer to (15, 20, 24 -26 33] for related discussion and surveys on 



some existing development. For more practical considerations, we shall work in a even 
more challenging setting by assuming very little a priori knowledge of the underlying 
scattering object, which might consist of multiple components, and both the number of 
the components and the physical property of each component are unknown in advance. 
This setting would make our current study even more highly non-trivial, but on the 
other hand would be of significant practical importance when there is little a priori in- 
formation available on the target object. Before we proceed to discuss more about our 
results, we next present the mathematical framework that we shall work within. 

Let O be a bounded C 2 domain in IR 3 which supports an inhomogeneous isotropic EM 
medium characterized by the electric permittivity e(x), magnetic permeability /J,(x), and 
conductivity a(x). Both e(x) and fi(x) are positive scalar functions and a(x) is a non- 
negative scalar function. It is assumed that Q e := IR 3 \f2 is connected and O e represents 
the uniformly homogeneous background space. Without loss of generality, we assume 
that f2 e is the vacuum, i.e., e(x) = fi(x) = 1 and a(x) = for x G VL e . But we would 
like to emphasize that all our subsequent studies could be straightforwardly extended to 
the case that f2 e is not necessarily the vacuum but uniformly homogeneous. The inverse 
problem that we shall consider in this paper is to recover the inhomogeneity (Q;e,n,a) 
by measuring the wave field far away from the inhomogeneity produced by sending a 
single detecting EM wave. In the physical situation, the inhomogeneity (f2; e, fj,, a) is 
referred to as a scatterer, and throughout we shall take the detecting/incident EM wave 
to be the time-harmonic plane wave, 

E\x)=pe iuJX - d \ H\x) = — V AE l (x), x G M 3 , (1.1) 

where u G IR+ denotes the frequency, A the exterior product, 6' G § 2 the impinging 
direction, and p G M 3 the polarization with p ■ 8' = 0. The EM wave scattering is 
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governed by the Maxwell equations, 



V A E - iuixH = 



a 



V AH + iu}(e + i-)E = 



111 

in 



E- 



E\ 



lim \x\ 



E+ = (E-E i )\ sle , 
(V A E + )(x) A t^t - iujE + {x) 



(1.2) 



0, 



where E and H are the respective electric and magnetic fields, E is the interior electric 



field, E + is known as the scattered electric field (cf. 15,28 ), and H~ and H + are 
their interior and scattered magnetic counterparts. Moreover, E + admits the following 
asymptotic development with the big O notaion 



E+( 



D iu)\x\ 



X 



-A 



1 



(1.3) 



x \x\ 



G S 2 



is 



uniformly in all directions x/\x\ as |x| —> +oo. A(8;9' ,p,u) with 6 
known as the electric far-field pattern. Here, we would like to distinguish two cases with 
< a < +oo and g = +oo. Following the terminologies in the literature, in the former 
case, (fi; e, p, g) is called an EM medium inclusion and the corresponding inverse problem 
is referred to as inverse medium scattering, whereas in the latter case as g —> +oo, the 
medium becomes perfectly conducting, and (0; e, fi, +oo) is called a perfect conductor or 
a perfectly electric conducting (PEC) obstacle, and the corresponding inverse problem 
is usually referred to as inverse obstacle scattering. For a perfectly conducting obstacle, 
both interior fields E~ and H~ := H\q would vanish inside f2, i.e., the EM fields cannot 
penetrate inside the object and are governed by 



V A E - iujH 

V A H + iuE 
vAE = Q 
E + = E-E i 



in O e , 
in f2 e , 
on dfl, 
in rh, 



(1.4) 



lim \x\ 

k |x|-s>+oo 



(VAE + )(i) A 



x 



icoE + (x) 



where v denotes the outward unit normal vector to d£l. The scattered field E + has 



the same asymptotic development as that in (1.3). Throughout the rest of the paper, 



we shall unify the notation of denoting a scatterer by (fi;e, n, g), which could be an 
inclusion by taking < g < +oo, or an obstacle by just taking g = +oo disregarding 
the parameters e and fi. 

The inverse scattering problem that we consider is to recover (fi; e, fj,, a) from the 
knowledge of A(9; 6',p, oS) or to recover fi from the knowledge of A(8; 8',p, oS) if g = +oo. 
We will make use of single far-field measurement, i.e., A{0; 9',p, uj) for all 6 G S 2 but 
fixed 6',p and u>. Furthermore, we shall require very little a priori knowledge of the 
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underlying scatterer, which could be composed of multiple components, and the number 
of the components is not required to be known in advance, and each component could be 
either a medium inclusion or an obstacle. Generically, in this extremely general setting, 
one cannot expect to recover all the details of the underlying scatterer by using only a 
single measurement. Instead, we would consider the recovery of locating the multiple 
components in the setting described above. Nevertheless, at this point, we would like 
to remark that our numerical experiments indicate that our proposed imaging scheme 
could also qualitatively reveal the supports/shapes of the scatterer components. 

Specifically, two single-measurement locating schemes would be proposed for two 
separate cases depending on the size of the target scatterer. The first scheme is for 
locating scatterers of small size compared to the wavelength of the detecting plane wave. 
The multiple scatterers could be extremely general with very little a priori knowledge 
required. Each scatterer component could be either an impenetrable perfectly conduct- 
ing obstacle or a penetrable inhomogeneous medium with an unknown content, and the 
number of the scatterer components is not required to be known in advance. The locat- 
ing scheme is based on a novel indicator function I s (z) for zeR 3 , If z happens to be 
the location point of a scatterer component, then z is a local maximum point for I s (z). 
Using the indicating behavior of I s (z), one could then locate all the scatterer compo- 
nents. In defining the indicator function, only inner-products of the electric far-field 
pattern and the vectorial spherical harmonics are involved. The indicating behavior is 
derived from the asymptotic expansion of the electric far- field pattern for small 'point- 
like' scatterers. The expansion is further based on the low frequency EM scattering 
asymptotics (i.e., Raleigh approximation); see [7}|9 p6|[27] . But for our imaging scheme, 
the expansion would be formulated in terms of the vectorial spherical harmonics instead 
of the polarizability tensors. It is also interesting to mention that our numerical experi- 
ments show that the proposed scheme also works effectively for locating 'partially-small' 
scatterer, namely, the scatterer is not 'point-like' but 'line-segment-like'. Furthermore, 
in addition to finding the locations of the scatterer components, the proposed scheme 
also shows some very promising feature in qualitatively imaging the supports/shapes of 
the unknown scatterers. 

Our second scheme is for locating multiple perfectly conducting obstacles whose sizes 
are comparable to the detecting EM wavelength. For this case, we would require that 
the shape of each component must be from a certain known admissible class. That 
is, there might be multiple obstacles with an unknown number of components, but the 
shape of each obstacle component must be from a certain class of reference obstacles 
which is known in advance. Nevertheless, there could be multiple different reference 
obstacles. Other than the assumptions made above, no further a priori information is 
required for the multiple unknown scatterers. The number of the unknown scatterer 
components could be larger than the number of the reference obstacles; that is, some of 
the unknown scatterers possess the same shape. Moreover, it is not necessary for all the 
reference obstacles to be presented in the unknown scatterer. The setting considered 
would be of significant practical interests, e.g., in radar and sonar imaging. The second 
imaging scheme is based on indicator functions, where I' £ N denotes the number of 
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the reference obstacles. Similar to the first imaging scheme, in calculating the indica- 
tor functions, only inner products are involved and no inversion will be required. The 
proposed method is very efficient and robust to measurement noise, and could also be 
extended to include inhomogeneous medium components if certain generic condition is 
satisfied. To our best knowledge, this is the first imaging scheme in the literature for lo- 
cating multiple unknown scatterers of regular size by a single EM far- field measurement. 
Rigorous mathematical justifications are provided for both schemes. In our theoretical 
analysis, we would impose the sparse distribution of the scatter components. However, 
we have conducted extensive numerical experiments, and numerical results show that 
even without the sparsity assumption, our proposed locating schemes still work very 
effectively. 

The rest of the paper is organized as follows. In Section 2, we develop the locating 
schemes, respectively, for the two separate cases with small scatterers and regular-sized 
scatterers. Section 3 is devoted to the proofs of Theorems 2.1 and 2.2 stated in Section 2 
to justify the indicating behaviors of the indicator functions for the two locating schemes 
proposed in Section 2. In Section 4, we present extensive numerical experiments to 
illustrate the effectiveness of the proposed methods. The work is concluded in Section [5] 
with some further discussions. 



2 The locating schemes 

In this section, we shall develop the two schemes for locating, respectively, multiple 
scatterers of small size and regular size. In order to ease the exposition, throughout 
the rest of the paper, we assume that oj ~ 1, and hence the size of a scatterer can be 
expressed in terms of its Euclidean diameter. 

2.1 Locating small scatterers 

We first introduce the class of small EM scatterers for our current study. Let I £ N and 
let Dj, 1 < j < I, be bounded C 2 domains in M 3 . It is assumed that all -Dj's are convex 
and contain the origin. For p 6 M + , we let pDj := {px\x £ Dj} and set 

0^ = Zj + pDj, zj eR 3 , 1 < j < I. 

Each £lj is referred to as a scatterer component and its content is given by Ej , pj and 
aj. It is assumed that £,• > 0,Pj > and Gj > are all constants, except the case that 
Oj = +oo. As remarked in Section [TJ if Uj is taken to be +oo, then Dj would be regarded 
as a perfectly conducting obstacle disregarding the parameters £j and pj. In the sequel, 
we shall reserve the letter I to denote the number of components of a scatterer defined 
by 

i i 
Q = fl, and s, p, cr) = [^J (Qj) £j, p,j, aj). (2-1) 
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The parameter p G M + represents the relative size of the scatterer (or, more precisely, 
each of its component). We now make the following qualitative assumptions, 



p < 1 and dist(zj, Zj>) > 1 for j / f, 1 < j,f < I. (2.2) 



The assumption (2.2) means that compared to the wavelength of the detecting/incident 
wave, the relative size of each scatterer component is small and if there are multiple 
components, they must be sparsely distributed. Our numerical experiments in Section 
4 could be more quantitative in this aspect, and it is numerically shown that if the size 
of the scatterer component is smaller than half a wavelength and the distance between 
two distinct components is bigger than half a wavelength, one could have a fine recon- 



struction by using our proposed scheme. In this sense, the qualitative assumption (2.2) 
is needed only for our subsequent theoretical analysis of the proposed locating scheme. 
Furthermore, we would like to emphasize that most of the other restrictive assumptions 
introduced above are also mainly for the purpose of the theoretical justification, and the 
numerical experiments in Section 4 will show that our proposed imaging scheme works 
in a much more general setting. Specifically, the regularity of each reference component 
Dj is not necessarily C 2 -smooth, and it could be a Lipschitz domain, and moreover, it is 
not necessarily convex and not even simply connected, as long as it contains the origin 
and pDj for a small p € M+ yields an appropriate domain of small size. Moreover, the 
content of each medium component is not necessarily constant, and it could be variable. 
Some of those restrictive assumptions could be relaxed from our theoretical justification 
in the following. However, it is our focus and emphasis on developing the locating scheme 
in the present work, and we would not appeal for a most general theoretical study. 

In the sequel, without loss of generality, we can assume that there exists I' G NU{0}, 
< 1' < I such that 

< cr j < +oo for 1 < j < I' and aj = +oo for V + 1 < j < I. 

Then, we set 

V I 

Q m ■= [J fij and n o := [J Qj (2.3) 

j=l j=l'+i 

to denote, respectively, the medium component and the obstacle component of a scat- 
terer. We emphasize that both / and /' are unknown in advance, and /' could be or 
/, corresponding to the case that O m = or Q D = 0. Moreover, the contents and the 
shapes of the two components are unknown in advance either and they are not even 
necessarily identical to each other, i.e., 

(Uj ; Ej , pj ,(Tj)^(n f ; e y ,p f , a r ) for j / f . 



Corresponding to a single EM detecting plane wave ( 1.1 ), the electromagnetic scattering 
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is governed by the following Maxwell system 



V A E - iuj | 1 + J2(fXj - l) X n/j H = 

V A H + + J2( £j - l)xn,-) - E ^ = 

v AE+ = -vAE { on dft , 

(V AE + )(x) A-^- -iojE + (x) = 0. 



in M 3 \fi c 



in K 3 \fi , 



(2.4) 



lim |x 

k \x\— >+oc 



It is known that there exists a unique pair of solutions (E+,H+) G # Zoc (curl; M 3 \ft ) A 
i?/ oc (curl; R 3 \f2o) to the Maxwell equations (2.4) (see [28]). In the following, we shall 
write 



A(9; n) := A(9; [j (fy ; ^ , ^ , ^ ) , 0', p, 



to denote the electric far-field pattern corresponding to the EM wave fields in (2.4). We 
shall also write for 1 < j < I, 

A{9;Q, j ) := A(9; (%; £j, aj), 0',p, ui) 

to denote the far-field pattern corresponding solely to the scatterer (flj) £j, /u, cr ? ) . Both 



15 



28 



A(9; 0) and A(9; Oj) are real analytic functions on the unit sphere S (cf. 

Next, we introduce the space of L 2 tangential fields on the unit sphere as follows, 



)• 



T Z (S Z ) := {a £ C | a G L 2 (§ 2 ) 3 , # • a = a.e. on § 2 }. 

Note that T 2 (§ 2 ) is a linear subspace of L 2 (§ 2 ) 3 , the space of vector L 2 -fields on the 
unit sphere S 2 . In addition, we recall the vectorial spherical harmonics 



1 



U™(9) := - , 

^Jn(n + 1) 

V™(9) :=9AU™{9) 



Grady™ (6) 



n G N, m 



-n 



,n, 



(2.5) 



which form a complete orthonormal system in T 2 (§> 2 ). In (2.5), Y™(0) with 9 £ § 2 , 
m = —n, . . . , n are the spherical harmonics of order n > 0, and Grad denotes the surface 
gradient operator on §. It is known that A(9; O) and A(9; JL-), 1 < j < I all belong to 
T 2 (S 2 ). We define 



K j :-- 



\\mnj)\\h>(s2) 

\\A(9;n)f T2(s2) ^ 



l<j<l, 



(2.6) 
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and 



h{z) :-- 



\\A(e ; n)f 



E 



T2(S2) m= _i,o,l 



> z u?(e) 



+ 



T 2 (S 2 ) 

A{9;n),e iuj{e '- e) - z V 1 m (e) 



(2.7) 



T 2 (§ 2 ) 



where z G M 3 , and (u, v)y2(§2) = J g2 u • v dsg. Clearly, is a real number, whereas 
is a real- valued function depending on the point z £ M 3 . We are now ready to 
present the first main result of this paper, whose proof will be postponed to the next 
section. 



Theorem 2.1. Let (fl;e, p, a) be given by (2.1) satisfying (2.2), and K 3 , l<j<l, and 
I s (z) be defined in (|2.6|) and (2.7), respectively. Set 



L = min distfz,-, Zy) 3> 1. 
i<M'<tjW 



Then 



K> - A;', + C? ( i + p 



i < i < z, 



(2i 



where Kq is a positive number independent of L and p. Moreover there exists an open 
neighborhood of Zj, neigh(zj), 1 < j < I, such that 



I s {z) <K J Q + o(^+p\ for ze neigh(zj), 
and moreover I s (z) achieves its maximum at z in neigh(zj), 

I B (zj) = K 3 +o(j+pY 



(2.9) 



(2.10) 



Remark 2.1. Clearly, I s {z) possesses the indicating behavior which could be used to 
identify the location point Zj's of the scatterer components Such behavior is more 

evident if one considers the case that f2 has only one component, i.e., 1 = 1. In the 
one-component case, one would have that 



I s (z) < 1 + O(p) for all z G M 3 \{.zi}, 



but 



I s ( Zl ) = 1 + O(p). 
That is, Z\ is a global maximum point for I(z). 



Using Theorem 2.1, we can formulate our first imaging scheme of locating multiple 
small scatterer components as follows, which shall be referred to as Scheme S in the rest 
of the paper. 

The first single-shot locating method: Scheme S 
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Step 1. For an unknown EM scatterer in (2.1), collect the far-field data by sending 
a single detecting EM plane wave specified by (1.1). 



Step 2. Select a sampling region with a mesh Th containing 0,. 
Step 3. For each point z 6 Th, calculate I s (z). 

Step 4. Locate all the significant local maxima of I s (z) on Th, which represent the 
locations of the scatterer components. 



2.2 Locating regular-sized scatterers 



In this section, we consider the locating of multiple scatterers of regular size. We shall 
present our scheme for the case that all the scatterer components are perfectly conducting 
obstacles. Later, we would remark how our method could be extended to include the 
inhomogeneous medium components. Let Mj C M 3 , 1 < j ' < I, be bounded simply 
connected C 2 domains that contain the origin. Let 



I i 
i=l j=i 



Zj e 



l < j < I. 



(2.11) 



Each Qj 



zj + Mj is a scatterer component, and it is assumed to be a perfectly 



conducting obstacle. Illuminated by a single incident EM plane wave (1.1), the EM 
scattering by the scatterer Q can be described by the following Maxwell system 



VAE- iujH = 
V A H + iujE = 
E+ = (E-E^ 
v A E + = —v A E 



in R 3 \Q, 
in R 3 \Q, 

on dQ, 



(2.12) 



lim \x\ 

Jx|-H-oo 



(V AE + ){x) A 



iu)E + (x) 



0. 



In the sequel, we write 



I i 
A(6; Q) := A(9; (J Qj) = A(9; (J ^ + Mj 

3=1 3=1 



(2.13) 



to denote the far-field pattern corresponding to the EM fields in (2.12). It is assumed 
that 

diam(Q j ) = diam(M,) ~ 1, 1 < j < I; (2.14) 

and 

L = min dist(^j, Zf) 3> 1. (2-15) 

l<3,3'<l,3j=3' 

That is, the size of the underlying scatterer components are comparable to the wave- 
length of the detecting EM plane wave. This is in sharp difference from our study in 
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Section 2.1 where the scatterer components are of small size compared to the detecting 



EM wavelength. The qualitative condition (2.15) states that the scatterer components 



are sparsely distributed, and we would like to emphasize again that this is mainly needed 
for our subsequent theoretical analysis of the proposed locating scheme. Our numerical 
examples in Section 4 shows that as long as the distance between different components 
are bigger than half a wavelength, the proposed scheme would yield a fine reconstruction. 
Furthermore, we introduce an admissible reference scatterer space 

? := {£,-}5U, (2-16) 

where each Sj C K 3 is a bounded simply connected C 2 domain that contains the origin 
and 

E^Sy/, for j^f, l<j,j'<l'. (2.17) 
For the present study, we assume that 



MjE^, j = 1,2,. ..,1. 



(2.18) 



That is, in the practical situation, the shapes of the underlying scatterers are known in 
advance in the sense that each component must be of a shape from a known admissible 
class. But we do not know the locations of those scatterer components and intend to 
recover them from a single wave detection. We would like to remark that it is not 
necessary to have 1 = 1'. It may have I > V , and in this case, there must be more than 
one component in Q who has the same shape from 5?\ and it may also have I < I 1 , and 
in this case, there are less scatterers presented than the known admissible scatterers. 
Next we introduce V indicator functions as follows, 



{A(9;Q),e i ^ e '- e > z A{e-^ k )) T ^ ) 

\\A(e-Y, k )\\l 2{s2) 



1,2, 



(2.19) 



where A(9; E&) denotes the far-field pattern corresponding to the perfectly conducting 
obstacle T* k , 1 < k < I'. In the following, we shall show that the I' indictor functions 



introduced in (2.19) can be used to locate the scatterer components Qj of Q. To that 



end, we shall make the following generic assumption that 

A(0; Efc) / A(9; E k ,), for k^k', 1 < k, k' < l'. 



(2.20) 



Assumption (2.20) is closely related to a longstanding problem in the inverse scattering 



theory (cf. [13,15,25]): whether or not can one uniquely determine an obstacle by a 
single wave measurement? That is, if two obstacles produce the same far-field data 
corresponding to a single incident plane wave, can one conclude that they must be 
the same? Though such uniqueness result is widespreadly believed to be true, it still 
remains open, particularly for an obstacle of regular size. The uniqueness is proved for 



obstacles of polyhedral type in 24,26 . Assumption (2.20) on the reference scatterers is 



of critical importance in our subsequently proposed locating scheme. Nonetheless, since 
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the admissible class is known, (2.20) can be verified in advance. Moreover, since 5? 
is known in advance, we can assume by reordering if necessary that 

||A(0;£ fe )|| T2(s2) > ||A(0;X fc+1 )|| r2(§2) , fc = l,2,.../-l. (2.21) 

That is, the sequence {||A(0; ^k)\\T 2 (§ 2 )Yk=i ^ s nonincreasing. Next, we present a key 
theorem, which is the basis of our subsequent locating scheme. 



Theorem 2.2. Let Q be given in (2.11), and the obstacle components are assumed to 
satisfy (2.15) and (2.18). The admissible reference scatterer space 5? is assumed to sat- 
isfy (2.20) and (2.21). Consider the indicator function 1} introduced in ( 2.19| ). Suppose 
for some jo G {1, 2, . . . , I}, Mj = Si, whereas Mj ^ Si for j £ {1, 2, ... , l}\{jo}. Then 
for each Zj, j = 1,2, ... ,1, there exists an open neighborhood of Zj, neigh(zj) , such that 



(i)- if j = jo, then 



ll(z) 



\lliz) - II < O 



L 



z E neigh(z 



'J0 7j 



(2.22) 



and moreover, Zj is a local minimum point for l}(z); 
(ii). if j ^ jo, then there exists eo > such that 



l}(z) := \ll(z) - 1| > e + O 



z G neigh(zj) 



(2.23) 



By using Theorem 2.2, our proposed locating scheme can be proceeded as follows, 
which shall be referred to as Scheme R in the rest of the paper. 

The second single-shot locating method: Scheme R 



Step 1. For an unknown EM scatterer Q in (2.11), collect the far-field data by sending 
the detecting EM plane wave specified by (1.1). 

Step 2. Select a sampling region with a mesh Th containing Q. 

Step 3. Collect in advance the far-field patterns associated with the admissible reference 
scatterer space 5? in (2.16), and reorder 5? if necessary to make it satisfy (2.21), 
and also verify the generic assumption (2.20). 

Step 4. Set ft = 1. 

Step 5. For each point z G Th, calculate I* (z) (or 1% (z) = \I^{z) — 1|). 

Step 6. Locate all those significant local maxima of 1% (z) such that I^{z) ~ 1 (or the 
minima of if (z) on Th such that I^{z) <C 1), where scatterer components of the 
form z + Sfc is located. 



Step 7. Trim all those z + Sfc found in Step 6 from Ti- 



ll 



Step 8. If Th = or k = V, then Stop; otherwise, set k = k + 1, and go to Step 5. 



It can be seen that our locating scheme R progresses in a recursive manner. For a 
scatterer Q with multiple components, one firstly locates the sub-components of shape 
Ei j which have the most prominent scattering effect among all the scatterer compo- 
nents. After locating all the sub-components of shape Ei, one can exclude those com- 
ponents from the searching region, and then repeats the same procedure to locate all 



the sub-components of shape E2, which, according to the ordering (2.21), have the most 



prominent scattering effect among all the scattering components that still remain in the 
searching region. Clearly, this procedure can be continued till one locates all the scat- 



ter components. Theorem 2.2 remains true for inhomogeneous medium scatterers if the 



generic condition (2.20) still holds and in that case, the locating Scheme R could be 



extended to include inhomogeneous medium components as well; see Remark 3.3 in the 
following for related discussions. 



3 Proofs of Theorems 2.1 and 2.2 



This section is devoted to the proofs of Theorems 2.1 and 2.2, which are the theoretical 
cores for our locating Schemes S and R, respectively. We first derive two key lemmas. 



Lemma 3.1. Let 



1 1 
r = |J Tj and (T;e, n,v) = \j(Tj;Ej, fij,<Tj), 
3=1 i=i 



(3.1) 



be a scatterer with multiple components, where each Tj is a bounded simply connected 
C 2 domain in M 3 . Assume that 



L = min dist(Tj, Tji) S> 1. 
i<j,j'<i,j¥=j' 



(3.2) 



Then we have 



(3.3) 



A{e-V) = Y J A{e;Y j ) + 0{L- x ), 
i=i 

where A(6;T) and A(8;Tj) denote the far-field patterns, respectively, corresponding to 
(T;e,n,a) and {Y y. ~ ; . n/j. 

Proof. We shall only consider a specific case with I = 2, and T± is an obstacle component 
while T2 is a medium component. Nevertheless, the general case can be proved following 
a completely similar manner. Moreover, we shall assume that oj is not an interior EM 



eigenvalue for T\ (see Remark 3.1). For this case, the EM scattering corresponding 



to r in (3.1) is governed by the Maxwell system (2.4) with f2 G replaced by T\ and 



Q, m replaced by (F^; £2, 1^2, 02)- We know that there exists a unique pair of solutions 
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E G (#/ oc (IR 3 \ri)) 3 and H G (F/ oc (M 3 \ri)) 3 to the Maxwell system (cf (28)). Next, we 
shall make use of the integral equation method to prove the lemma. To that end, we let 

9(x,y):=—- r, x,y£R s , x^y, 

47r \x — y\ 

which is the fundamental solution to the differential operator —A — uJ 2 . Let R G M+ be 
sufficiently large such that ri U T2 <s -Br. 

By the Stratton-Chu formula [15], we have 



and 



E(x) =E\x) + V x A / u(y) A E(y)$(x, y) ds y 

JdT! 

- —V x A V x A / A H(y)<f>(x, y) ds y 

+ V* A / (Vy A E(y) - iujH{y))$(x, y) dy 
Jv 2 

- V x f VyE(y)$(x,y) dy 

+ iu (V y AH(y) + iujE(y))<S>(x,y)dy, x G B R \T U 
Jr 2 

H{x) =H i (x) + V x A / v(y) A H(y)$(x, y) ds y 
JdT 1 

+ — V x A V x A f v(y) A E{y)®{x, y) ds y 

+ V X A (VyAH{y)+iojE(y))^(x,y)dy 
Jr 2 

- V* / VyH{y)${x,y) dy 

Jr 2 

-iu (V y AE(y)-iuH(y))^{x,y)dy, x E B R \Tx. 



r 2 



(3.4) 



(3.5) 



By using the Maxwell equations (2.4), and (3.4)-(3.5), together with the use of the 



mapping properties of the layer potential operators (cf. [151128]) , we have the following 
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system of integral equations 

E(x) =E\x) - —V x AV,A / v(y) A H{y)${x, y) ds y 

+ V X A iu((i 2 - l)H(y)$(x, y) dy 
Jt 2 

+ ioj(l-e 2 -i^E(y)$(x,y) dy, x£T 2 , (3.6) 

H(x) =H\x) + V x A / u(y) A H(y)$(x, y) ds y 
JdT! 

+ V X A iu)(l-e 2 - i—) E(y)$(x, y) dy 
Jt 2 v u J 

-[ ioj(ji 2 -l)H(y)<S>(x,y) dy, x G T 2 , (3.7) 
Jt 2 

v(x) A H(x) =2v(x) A H\x) + 2u(x) AV,A / v(y) A H(y)$(x, y) ds y 

JdT! 

+ v(x) A V x A / 2iu (l - e 2 - i—) E{y)<$>{x, y) dy 
Jt 2 v 00 ' 

- v{x) A / 2iu(fi 2 - l)H(y)Q(x, y) dy, x G dT±. (3.8) 
Jt 2 

Next, we introduce the following volume integral operators 

(CE)(x) = [ iuj (l - e 2 - i—) E{y)$>{x, y) dy, x G T 2 , 
Jt 2 v w 7 

(£'E)(x)=V x A(£E)(x), xGT 2 , 



(K,H)(x) = / iuiji* ~ l)H(y)$(x,y) dy, x G T 2 , 
Jt 2 

(K'H)(x) =V X A (ICH)(x), x G T 2 , 
and the boundary integral operator 

^P{y A H))(x) = 2v{x) A V x A / AH(y)$(x,y) ds y , a; € STi. (3.10) 

Two important function spaces are recalled as follows: 



T#- 1/2 (dr!) := | a G (#~ 1/2 (dri)) 3 ; i/Aa = for a.e. s e ^j} , 

TH mJ 2 ( dr i) : ={ a G ^-^(anj^ivCa) g rr 1 / 2 ^!)} . 



(3.11) 
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Furthermore, we let 



bi(x) :=--V,AV,A / u(y) A H (y)$(s, y) ds y , x E T 2 , 

b 2 (x) :=V X A v(y) AH{y)$(x,y) ds y , x E T 2 , 

b 3 (x) :=i/(a?) AV,A / 2iw f 1 - e 2 - i— ) E(y)${x, y) dy 

- i/(ar) A / 2iw( / u 2 - l)if(y)$(x, y) dy, x E IY 
Jr 2 



T 2 

Since L = dist(ri,r2) 1, one readily verifies that 

C, 



) 'llL 2 (r 2 )3 



< — hvAi? 



rffDi L v /2 ( 9r 



,1 = 1,2, 



and 



C 



H b 3ll TH -i/2 (ri) < j (ll^lli 2 (r 2 ) 3 + \\ h \\l^(T2) 3 ) > 
where C is a positive constant depending only on ri,r 2 and u. 



(3.12) 



(3.13) 



Next, by using the integral operators introduced in (3.9) and (3.10), the integral 



equations (3.6)-(3.7) can be formulated as 
A :: 



I-C K 1 



E(x)\ fhi(x)\ _ (E l (x 



W(x) ' 



H(x)J \b 2 (x) 
In a similar manner, ( |3.8[ ) can be formulated as 

(/ - V)(y A H)(x) - b 3 (s) = 2u{x) AH\x), x E r x . 



x E T 2 , (3.14) 



(3.15) 



Referring to [15 28 , we know that both A : L 2 (r 2 ) 3 A L 2 (r 2 ) 3 — > L 2 (r 2 ) 3 A L 2 (r 2 ) 3 and 

— \/2 —1/2 

I — V : TH Di ^ (dTi) — > TH Di ^ (dFi) are Fredholm operators of index 0, and thus they 
are invertible. Using this fact and ( |3.12[ )— (3.13 ), one obtains from (3.14) that 



E(x) 
H(x) 



-i 



E\x) 
HHx 



+ 0(L 



E(x) 
H{x) i 



+ D{L 



x e r 2 



(3.16) 



and 



(vAH)(x) = (I-Vy 1 (2uAH i )(x)+0{L- 1 ) := (y AH)(x) +0{L~ 1 ), x E dT 1} (3.17) 

where (E,H) are actually the EM fields corresponding to (r 2 ; e 2 , /i 2 , cj 2 ), and -ff is the 
magnetic field corresponding to T\. Finally, by using (3.16) and (3.17), and the integral 

□ 



representation (3.6), one readily has (3.3). 



The proof is complete. 
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Remark 3.1. If a; is an interior eigenvalue for Ti, one can overcome the problem by using 
the combined electric and magnetic dipole operators technique, and we refer to [l5| Chap. 
6] for more details. 



Lemma 3.2. Let (Qi; e%, p,\, a\) be one component of£l described in Section 2.1 . Then 
we have 



A(0-Sl x ) =e 



iu(6'-0)-z 1 



A(9- P D l 

m=-l,0,l 



(3.18) 



where U™ and V™ are the vectorial spherical harmonics introduced in fl2.5|) ; and axm 
and bi >rn (m = —1,0,1), are constants depending only on (D\; e±, pi, a\) and p,6', but 
independent of cop. 

Proof. We first consider the case that a\ = +oo, namely, J7i is a perfectly conducting 
obstacle. It is directly verified that 



A(0; fi x ) = A(9; z\ + pD x ) = e l ^ 6 '- 6 > z ^A(6; pD{). 
The EM scattering corresponding to the obstacle pD\ is described by 



(3.19) 



V A E p - iujH p = 

V A Hp + iujEp = 
E+ = (E p - E? 



vAE+ = 
lim \x\ 

k |z|->+oo 



-v A E l 



(VAEt)(i) A — -iuE+{x) 



in R d \pDi, 
in R 3 \p~Th, 

on d(pD{), 
= 0. 



(3.20) 



Set 



E{x) := E p (px), H(x) := H p (px), E i {x) = E i {px) for x € M 3 \L>i 
It is verified directly that 

'\7AE-iujpH = in R 3 \L\, 



(3.21) 



V A H + iwp-E = 
E + = {E - E l )\^^, 
v A E+ = -v A E i 



in M 3 \D 1 



on <9-Di, 



(3.22) 



lim \x\ 

k |a;|— >+oo 



(V A £+)(x) A —r - iw£7 + (x) 



0. 



Then by the low frequency asymptotics in [lfjj, one has 



£ + (:e) 



iup\x 



A(8) + O 



(3.23) 
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with 



A(6) = ^^[9 A (9 A a) - 9 A b] + 0{{ujpf), 



(3.24) 



where a and b are two constant vectors, representing the electric and magnetic dipole 
moments, and they depend only on p, 9' and D\, but independent of to and p. By (3.21 ), 
one can readily see that 



A{9) = —A{9) = JL ("M! [# A (0 A a ) _ e A b l + 0((up) 4 ). 
iuj iuj Air 



(3.25) 



Finally, by using the fact that U™ and V™ form a complete orthonormal system in 
T 2 (§ 2 ), it is straightforward to show that there exist a\^ m and &i jm , m = —1,0, 1 such 
that 

#A(#Aa)-0Ab= Yl ^,mur(0) + h, m vr(9), 

m=-l ,0,1 



which together with (3.24) and (3.19) implies (3.18). 

For the case when the underlying small scatterer is an inhomogeneous medium, by us- 
ing a completely same scaling argument, together with the corresponding low frequency 
EM asymptotics in 16 , one can prove (3.18). 

The proof is completed. □ 

Remark 3.2. The low frequency EM asymptotics and the asymptotic expansions of EM 
fields due a small inclusion could also be found in [7[|9||27|, where polarizability tensors 
are always involved. For the present study, we need asymptotic expansions in terms of 
the vectorial spherical harmonics. 

After technical preparations above, we are now in a position to show the proof of 
Theorem EH 



Proof of Theorem 2.1. First, by Lemmas 3.1 and 3.2, we have 



[ /i 
A(9;n) = ^A(9;n l ) + 0(- 

.7=1 ^ 



I 



Y^ e M0'-0)-Zi [(ujpfA j (9) + O((L0pf)] +0 



iu){9'~6)-z 

i=i 
i 

i 



m =-l,0,l 



1 

L 



(3.26) 



+ o 



L 



where ai m and b^ m , m = —1,0, 1, are constants dependent only on (Dj;ej,pj,aj) and 



p, 9'. In (3.26), we have introduced 



A'"(0):= £ {a^ m Ur(9) + b^ m Vr(9) 
m =_i,o,i 
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Next, without loss of generality, we only consider the indicating behaviors of I s (z) 
in a small open neighborhood of z±, i.e., z S neigh{z\). Clearly, we have 



UJ\Z 



■j — z\ > ujL S> 1 for z E neigh(zi) and j = 2, 3, . . . , I. 



(3.27) 



Hence, by the Riemann-Lebesgue lemma and (3.26), we have 
(A(9;n),e^ e '- e ^Ur(0)) T 2 m 



(3.28) 



<M 3 ( a fi + o[ 1 + P 



m 



■1,0,1. 



In (3.28), we have employed the orthogonality of vector spherical harmonics and also 
the Cauchy-Schwartz inequality. Moreover, by the Cauchy-Schwartz inequality, we know 
the strict inequality would hold in the last inequality of (3.28) if z ^ z\, and only when 
z = z\, the equality would hold. In a completely similar manner, one can show that 

=(upf{e i ^ e '- d ^A\e), e^'-^VTWW*) + O ( (upf{\ + P) 



L 



(3.29) 



m 



-1,0,1, 



where strict inequality would hold for the last relation if z 7^ Zi, and only when z = Zi, 
the equality would hold. 

Hence, by (3.28) and (3.29), we have 



E 

m=-l,0,l 



(A(0:n), ( ^°'- >*Ur(0) ri 



+ 



A(e;n),e iw( - e '- e > z vr(d) 



T 2 (§ 2 ) 



(3.30) 



<(-p) 6 ( E \%1\ 2 + \^+°(t+p)) 



where the equality would hold only when z = z\. On the other hand, by using (3.26), it 
is straightforward to show that 



\\A(9;n)r THS2) = (upfY,[ E wU 2 +Kl\ 2 + [i+p 



j=l \m=-l,0,l 



1 



(3.31) 
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By (3.30) and (3.31), we see that for z 6 neigh{z\) 
1 



E 



T 2 (S 2 ) m =-l,0,l 



A(0;fi),e 



lOJItf —0)-Z T TTfl 



U^{6) 



+ 



T 2 (§ 2 ) 

A(e;Q),e iu,i9 '~ 9) - z v; n (e) 



T 2 (§ 2 ) 



1 



E l a l,m| 2 + l^l,m| 2 

< -^5==* +0 



/ 1 



E E K' 



+ p 



j'=l m=— 1 

where only when z = z\, the equality would hold in the last relation. Set 



Kl,- 



E l«£| 2 + lCl 2 

m=— 1 

E E l«Sll 2 + l^il 2 

j'=l m=— 1 



Using (3.18) and (3.31), it is readily seen that 

\\A(9-M\\^ (Sl) _ 



K 1 



\\A(9;mh(S2) 
Thus the proof is completed. 



Kl + 0[ 1 + p 



Next, we show the proof of Theorem 2.2 



Proof of Theorem 2.2. First, by Lemma 3.1 we have 

l l 



A(0; Q) = A(9; \J Qj) = £ A(9; Qj) + O ( i 

j=l j=l V 

.7=1 .7=1 V J 



3=1 



(3.32) 



□ 



(3.33) 



Let us consider the indicator function l}{z) in (2.19). Without loss of generality, we 
assume that M\ = Si and Mj ^ Si for j = 2,...,/. Let z £ neigh{z\). By (2.15), 
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(3.33) and the Riemannian-Lebesgue lemma, we have 



A(#; s i))t 2 (s 2 ) 



+ o 



L 



(3.34) 



where in the last relation, the equality would hold only when z 
sufficiently small neighborhood, neigh(z{), 

1 



z\. Hence, for a 



\I r (z) - 1| < O 



L 



for z G neigh(zi). 



Next, we let z G neigh(z2), with neigh{z-i) sufficiently small. Then, again by (2.15), 
(3.33) and the Riemannian-Lebesgue lemma, we have 

(A(e;n),e^ e '- e >^A(e^ 1 )) T2i§2) 



(§2) 



+ o 



(3.35) 



<||A(0;M 2 )|| T2(S2) ||A(^S 1 )|| T2(S2) +O 



where in the last relation, strict inequality holds in light of the assumption (2.20) and 
the fact M2 7^ Si. Hence, by (3.35) and (2.21), it is readily shown that 

1 



^•(22) < -n-r77r-^m — — < 1 + O 



Hence, there exists some eo G M + such that 

1 

z 



(3.36) 



/^(z) = \ll(z) - 1| > e + O 



for z G neigh(z2) 



In a completely similar manner, one can show the indicating behaviors of l}{z) for 
z G neigh(zj), j = 3, ... ,1. 

The proof is completed. □ 



Remark 3.3. Through our proof of Theorem |2.2[ one can see that our locating Scheme 
R could be extended to a more general setting by including inhomogeneous medium 
components as follows. Suppose in y, some reference scatterer, say, (Ej-; £j, [ij, aj), 
1 < j < /' < I, are inhomogeneous media, with the EM parameters £j,^j and dj known 
as well. Let Q be a scatterer with multiple components, such that each component is a 
translation of some reference scatterer. Then our locating Scheme R could be extended 
to this more general setting provided the reference scatterer space ,y satisfies the generic 
assumption (2.20). However, to our best knowledge, there is no such uniqueness result 
in the literature. On the other hand, as we remarked earlier that y is given in advance, 
one coud verify (2.20) in advance, and if it is satisfied, our locating scheme could apply. 
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4 Numerical experiments and discussions 



In this section, we have carried out a series of numerical experiments for different bench- 
mark problems to test the performance of our proposed locating Schemes S and R. 
The results achieved are consistent with our theoretical predictions in Sections [2] and [3] 
in a sound manner. Besides, the numerical results reveal some very promising features 
of the imaging schemes that were not covered in our theoretical analysis. 

We first briefly describe our experimental settings. Let e\ = (1, 0, 0) T , e 2 = 
(0, 1, 0) T and = (0, 0, 1) T be the three canonical Cartesian bases. The single de- 
tecting/incident wave we shall employ for our numerical examples is the plane wave 



specified by (1.1) with the polarization p = and impinging direction 9 = e\. More- 
over, We shall take the unitary wavelength A = 1, namely the frequency u = 2n. In 
all the examples, the electric far-field pattern A is observed at 590 Lebedev quadrature 



points distributed on the unit sphere S (cf. 22 and references therein). The exact 



far-field data A{9) are corrupted point-wise by the formula 

A s {6) = A{9) + 5Cimax\A(6)\ exp(i2< 2 ) , (4.1) 
o 

where 5 refers to the relative noise level, and both (j and £2 follow the uniform distri- 
bution ranging from —1 to 1. The scattered electromagnetic fields are synthesized using 
the quadratic edge element discretization in the spherical domain centered at the origin 
with radius 4A enclosed by a spherical PML layer of width A to damp the reflection. 
Local adaptive refinement techniques within the inhomogeneous scatterer are adopted 
to enhance the resolution of the scattered field. The far-field data are approximated by 
the integral equation representation 29, p. 181, Theorem 3.1] on the sphere centered 
at the origin with radius 3.5 using the numerical quadrature. We refine the mesh suc- 
cessively till the relative maximum error of successive groups of far-field data is below 
0.1%. The far-field patterns on the finest mesh are used as the exact data. The values 
of the indicator functions have been normalized between and 1 to highlight the po- 
sitions identified. The sampling domain is fixed to be T = [— 2A, 2A] 3 , which is then 
divided into small cubes of equal width h = 0.01 A yielding the sampling mesh Th- The 
orthogonal slices of the contours of the indicator function values will be displayed as an 
estimate to the profiles of the unknown scatterers. In the sequel, for brevity, we shall 
refer to the newly proposed single-shot locating methods of Schemes S and R as the 
SSM(s) and SSM(r), respectively. 

The scatterers under concern include a cube, a ball of different radii, and revolving 
solids from a kite and a peanut parameterized in the x — y plane as follows [23] : 

(x(t),y(t),0) := (cos t + 0.65 cos It- 0.65, 1.5 sin i,0), < t < 2vr , 
(x(t),y(t),0) := y/3 cos 2 1 + l(cost, smt, 0), < t < 2vr . 

In the following, three groups of experiments shall be conducted. The first group of 
experiments is on locating point-like small-sized scatterers in various scenarios by the 
SSM(s), and the second group of experiments is on testing the SSM(r) for locating 
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Figure 1: Reconstruction results for Example 1: (Left) Exact far- field data; (Right) 
Noisy far- field data with 5 = 20%. 



regular-sized scatterers. In the third group of experiments, we shall test the performance 
of SSM(s) on imaging 'partially-small' line-segment-like scatterers. 

4.1 The SSM(s) for point-like scatterers 

Example 1. In this example, we consider a cube scatterer of length 0.02A located at 
the origin, with the EM parameters given by e = 4, \x = 1 and a = 0. The orthogonal 



slices of the contours of the indicator function I s (z) in (2.7) for the SSM(s) are given 
in Fig. [I] It can be readily seen that the SSM(s) can locate the small scatterer in a 
very accurate and stable manner. Indeed, even 20% random noise is attached to the 
measurement data, the SSM(s) still yields a very robust and accurate locating. 

Example 2. In this example, we consider a scatterer consisting of a ball medium com- 
ponent of radius 0.2A positioned at (1.5A, 1.5A, 0) and a ball PEC obstacle of radius 
0.2A positioned at (— 1.5A, — 1.5A, 0). The EM parameters of the first scatterer compo- 
nent are taken to the same as the one in Example 1. For this example, the orthogonal 
slices of the contours of the indicator function I s (z) for the SSM(s) are shown in Fig. [2] 
The SSM(s) yields a very accurate identification of the location of both scatterers even 
if the measurement data is significantly perturbed to a high-level noise. This exam- 
ple demonstrates that SSM(s) can locate the multiple scatterer components without 
knowing the physical property of each component in advance. Moreover, we would like 
to note a promising feature of the SSM(s): if one chooses a cut-off value to be 0.7 to 
separate the region where I s (z) is bigger than the cut-off value, then the rough profiles 
of the scatterer components would appear. Hence, in addition to find the locations of 
the scatterer components, our proposed scheme could also be able to qualitatively image 
the supports/shapes of the unknown scatterers. 

Example 3. In this example, three scatterers are presented at (—A, —A, A), (A, A, 0), 
(A, A, —A), respectively. The first two scatterers are balls of radius 0.2A, with the EM 
parameters given by e\ = 4 + sinxi, pL\ = 1, a\ = and £2 = 4 + cosx2, [ii = 1, 02 = 0, 
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Figure 2: Reconstruction results for Example 2: (Left) Exact far-field data; (Right) 
Noisy far- field data with 5 = 20%. 





Figure 3: Reconstruction results for Example 3: (Left) Exact far-field data; (Right) 
Noisy far-field data with S = 20%. 



respectively. Here, we have made use of x = (x\, x^-, x%) to denote a point in R 3 . The 
third one is a non-convex kite-shaped revolving scatterer scaled by a relative size p = 0.2 
with the EM parameters given by £3 = 4 + x$, //3 = 1 and 03 = 0. The numerical 
reconstruction results are shown in Fig. [3j It can be seen from Fig. [3] that the SSM(s) 
is capable of locating multiple scatterer components with variable contents. Moreover, it 
verifies that the small scatterers are not necessarily convex as assumed in our theoretical 



justification (cf. the remark below Eq. (2.2)). 

Example 4. We consider two ball scatterers of radius 0.2A, located at (— 0.45A, 0, 0) 
and (0.45A, 0, 0), respectively, with the same EM parameters e = 4, /i = 1 and a = 0. 
We shall investigate the lower distance limit between the underlying separate scatterer 
components for the SSM(s). The results are shown in Fig. |4j It can be seen in this case, 
namely the distance between the two components is of a half wavelength, the SSM(s) 
can locate both scatterer components and separate them well. If we further reduce the 
distance between the two components (less than a half wavelength), the SSM(s) can no 
longer separate the two scatterer components, though it can still roughly locate them. 
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Figure 4: Example 4: Reconstruction results for Example 4: (Left) Exact far-field data; 
(Right) Noisy far-field data with 5 = 20%. 

4.2 SSM(r) for regular-sized scatterers 

In this subsection, we shall consider two examples to demonstrate the capability and 
effectiveness of the SSM(r) for locating regular-sized scatterers. 

Example 5. In this example, we consider a revolving peanut and a revolving kite, both 
revolving solids along the x-axis. The gravitational center of the peanut is chosen at 
(0, 0, 1.5A), and the center of the kite is anchored at (0, 0, — 1.5A). Both scatterers are 
taken to be PEC obstacles. 

The admissible reference class is chosen to be composed of a unit ball, a revolving 
kite and a revolving peanut. All those admissible reference obstacles are centered at the 
origin and are PEC obstacles. Following the SSM(r) algorithm, their norms of far-field 
data associated with the reference obstacles in the admissible space are sorted in the 
descending order with the first one being the kite, the second one being the peanut 
and the third one being the unit ball. Using those a priori admissible known far-field 
data, we implement the SSM(r) method. The orthogonal contour slices with a certain 
transparency are shown in Fig. [5] for better visualization. For both the noise-free and 
noisy far-field data, the SSM(r) can successfully determine the location of the kite 
through the first indicator function; see the dark red part in the center of the kite in the 
left figures. Once the kite is determined and its surrounding subregion trimmed from the 
sampling domain, we find that the second indicator function plot shows that the center 
of the peanut can be identified by continuing the SSM(r), see the right figures. After 
the peanut is determined and trimmed from the sampling domain, if one continues the 
SSM(r) by testing data of the reference ball obstacle, no significant peak value would 
be found for the reference unit ball obstacle. Our numerical results are nicely consistent 
with our theoretical predication. 

Example 6. In this example, we adopt the same setting as that in Example 5, but with 
all the scatterer components set to be inhomogeneous media, with the EM parameters 
e = 4, [x = 1 and a = 0. It can be verified that the far-field data for the three reference 



24 






Figure 5: Example 5: Reconstruction results using (top) exact far-field data, (bottom) 
noisy far-field data with 5 = 5%. 



scatterers are distinct from each other, hence the generic condition (2.20) is satisfied and 



thus the SSM(r) applies; see Remark 3.3 Different from the obstacle case, the order 



of the norms of the far-field data of medium components in the admissible reference 
space in Example 6 is as follows, the peanut comes first, and then the kite and finally 
is the unit ball. The numerical results are given in Fig. [6j from which we can see that 
the SSM(r) can also successfully determine the locations of each medium component 
successively. 



4.3 SSM(s) for line-segment-like scatterers 

In this section, we present some numerical examples to show some very interesting and 
promising features of the SSM(s) that were not covered in our theoretical analysis. 
Specifically, we shall test the performance of the SSM(s) in identifying 'partially small' 
line-segment-like scatterers. 

Example 7. This example considers a slender cylinder scatterer with base point an- 
chored at (— 0.5A, 0, 0), radius 0.1A, height 2A and pointing to the positive x-axis. The 
EM parameters inside the slender cylinder are chosen to be e = 4, /i = 1 and a = 0. The 
numerical results are given in Fig. [7J In the noise-free case, the SSM(s) can determine 
successfully the location and even the length for the slender cylinder scatterer. The 
identified geometry is roughly a red long bar with the correct length. Furthermore, the 
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Figure 6: Example 6: Reconstruction results using (top) exact far-field data, (bottom) 
noisy far- field data with 5 = 5%. 

SSM(s) also performs well for the even more challenging case with large random noise, 
up to 5 = 20%, attached to the far-field data. 

Example 8. In this example, we consider an L-shaped scatterer composed of two 
slender cylinders as given in Example 7, except that both base points are moved to 
(—A, —A, 0) and pointing to the positive x- and y-axes, respectively. The EM parameters 
inside the scatterer are chosen to be e = 4, fj, = 1 and a = 0. For this example, it 
turns out that we need multiple probing wave measurements in order to have a fine 
reconstruction. Specifically, we shall take three far-field measurements corresponding to 
(Pi>#'i) = (e3,ei), (p2,G' 2 ) = (ei,e2) and (p3,0' 3 ) = (e2,es). The indicator function that 
we shall use for the identification is given by taking the maximum value of the three 
separate indicator functions corresponding to the three far- field measurements. That is 

I s (z) = max!<K3 4 l \z), z€T, 

where iP (z) is the indicator function calculated by using the far-field data generated by 
the plane wave with (pi, 9[),l = 1, 2, 3. It can be seen from Fig. |8]that an L-shaped dark 
red bar is identified from the composite indicator function. The reconstruction result of 
the SSM(s) is even robust to noise up to 20% as well. 
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Figure 7: Example 7: Reconstruction results using (left) exact far-field data, (right) 
noisy far-field data with 5 = 20%. 




5 Concluding remarks 



In this paper, two inverse scattering schemes, SSM(s) and SSM(r), are proposed for 
locating multiple electromagnetic scatterers by a single electric far-field measurement. 
The locating schemes could work in an extremely general setting: the underlying scat- 
terer might include, at the same time, obstacle components and inhomogeneous medium 
components; the number of the scatterer components and the physical property of each 
component are not required to be known in advance. The first scheme SSM(s) is for 
locating scatterers of small size compared to the detecting EM wavelength. For this 
scheme, the content of each medium component is not required to known in advance 
either. The second scheme SSM(r) is for locating scatterers of regular size compared 
to the detecting EM wavelength. For this scheme, certain a priori information would 
be required of each scatterer component. Specifically, if the component is an obstacle, 
then its shape must be from a certain known admissible reference scatterer class; and 
if the component is an inhomogeneous medium, then its support and content must also 
be from a certain known admissible reference scatterer class, and moreover, a certain 
generic condition must be satisfied. Nevertheless, the scheme SSM(r) could also work 
in a very general setting. The reference class may consist of multiple different refer- 
ence scatterers, and some reference scatterer may not be presented as a component in 
the unknown scatterer, and some may be presented as components for more than one 
time. The setting considered would be of significant interests, e.g., in radar and sonar 
imaging. The locating schemes are based on some novel indicator functions, whose in- 
dicating behaviors could be used for identification in many applications. In calculating 
the indicator functions, no inversions would be involved, so the proposed methods are 
very efficient and robust to noise. Rigorous mathematical justifications are provided for 
both schemes. Extensive numerical experiments are conducted to illustrate the effec- 
tiveness and robustness of the proposed imaging schemes in various practical scenarios. 
The numerical results match our theoretical predications in a very sound manner. Fur- 
thermore, the numerical results also reveal some very interesting and promising features 
of the proposed schemes. The scheme SSM(s) is also capable of qualitatively imaging 
the supports/shapes of the unknown scatterers in addition to locating them. Second, it 
is also capable of qualitatively identifying 'partially small' line-segment-like scatterers. 
Both of these topics are worth further investigation. 
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